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ABSTRACT 
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Context. Observationally, spectra of brown dwarfs indicate the presence of dust in their atmospheres while theoretically it is not clear 
what prevents the dust from settling and disappearing from the regions of spectrum formation. Consequently, standard models have 
to rely on ad hoc assumptions about the mechanism that keeps dust grains aloft in the atmosphere. 

Aims. We apply hydrodynamical simulations to develop an improved physical understanding of the mixing properties of macroscopic 
flows in M dwarf and brown dwarf atmospheres, in particular of the influence of the underlying convection zone. 
Methods. We performed two-dimensional radiation hydrodynamics simulations including a description of dust grain formation and 
transport with the C05B0LD code. The simulations cover the very top of the convection zone and the photosphere including the dust 
layers for a sequence of effective temperatures between 900 K and 2800 K, all with log g=5 assuming solar chemical composition. 
Results. Convective overshoot occurs in the form of exponentially declining velocities with small scale heights, so that it affects 
only the region immediately above the almost adiabatic convective layers. From there on, mixing is provided by gravity waves that 
are strong enough to maintain thin dust clouds in the hotter models. With decreasing effective temperature, the amplitudes of the 
waves become smaller but the clouds become thicker and develop internal convective flows that are more efficient in transporting and 
mixing material than gravity waves. The presence of clouds often leads to a highly structured appearance of the stellar surface on 
short temporal and small spatial scales (presently inaccessible to observations). 

Conclusions. We identify convectively excited gravity waves as an essential mixing process in M dwarf and brown dwarf atmo- 
spheres. Under conditions of strong cloud formation, dust convection is the dominant self-sustaining mixing component. 
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1 . Introduction reaching into ever cooler and lower mass regimes (Teff > 600 K, 

_ i M > 5 Mjup). 

Brown dwarfs form like stars and evolve as they cool from 
O stellar-like properties -M spectral type characterized by molec- Atmospheric temperatures are sufficiently low (T,.,, < 
5 ular hydrogen and water vapor formation, chromospheric ac- ^^OOK, Teff < 2800 K) for dust particle formation to occur in 
: tivity, flares, and magnetic spots - to planet-like properties - ^^^^-^yP^ M dwarfs. These grains should sink under the influ- 
> . T specti-al type characterized by methane absorption, electron- ^"'^e Suavity (g * 10^ cms-^) into deeper layers and van- 
■C^ • degenerate core, and maser emission. With a fully convective from the atmosphere, clearing it from condensable mate- 
<N ■ interior reaching up to the atmosphere, and a neutral atmosphere "^1- However, their near-infrared specti-a can only be repro- 
^ offering little interaction with magnetic field lines, they retain ^uced when accounting for a strong greenhouse eff^ect (also 
■ ■ ■ larger rotational velocities (> 30 kms-i, i.e., P < 4 hrs com- '^^"^'^ ^ blanketing eff^ect m steflar physics) m the visible lay- 
pared to 11 hrs for Jupiter). This efficiency and the unusually er^Tsun et al. 1996; Alexander et al. 1997; Ruiz et al. 1991 
large extent of the convection zone into the atmosphere (up to iLeggett et alj|1998t |2001). Classical static model atmospheres 
optical depths of 10 ^ for M dwarfs) assigns an important role as ^ave to rely on ad hoc assumptions about the mechanism that 
a cooling and contraction evolution regulator to the atmospheres ^eeps dust from settling, or that brings fresh materid towa rd the 
(lBaraffe etal.irT995h . Understanding the atmospheric properties ^'^i*'^^ grains to form (Helling et^|2008a). The 
has therefore implications for the mass determination of these ^^^'^^s °f ^^^^^ formation on the atmospheres of late-type dwarfs 
objects. Some 700 brown dwarfs have been founcO in the solar have been explored by modeling dust formation in chemical 
neighborhood and in star-forming regions since the eai-ly 90's equihbrium with Ae gas phase using diverse prescriptions of the 
^ * * ^ cloud thickness dAfl ardetal. 2001; Tsuu 2002; B urrows et al.1 
|2006'; 'A ckerman & Ma rlev 2001). It has been found that the 



Send offprint requests to: B.FieytSLg PHOENIX Dusty models (Allardetal. 2001) reproduce the in- 

1 Photometry spectroscopy and astrometry of M, L, and T dwarfs: frared emission of late-type M to mid-L dwarfs (i.e., 170 K < 
[ht^:/ /spideripac.caltech.edu/stafl7davy/ARCHI VE/index.shtml| T^g < 2500 K) dLeggett et al.ll998ll200lllRuiz et al.ll997i) . This 
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indicates that dust forms close to equihbrium in the infrared- 
line-forming region of these atmospheres (t w 10"^). 

On the other hand, late-type L and T dwarfs (reff < 1400 K) 
are less affected by photospheric greenhouse effects, but do show 
evidence - in terms of higher CO and lower NH3 and CH4 
absorption - of the dynamical upxyelUng of N 2 and CO gas 
("Saumo n etaTll2006h . ICushing et all ('2008) and IStephens et aP 
(2009) fitted a sequence of the red optical to mid-infrared spec- 
tra of_earl}ML_jo_jnid£rdwffl'fs us ing the model atmos pheres 
of lAckerman & MarlevI ^Wl) and' Saumon et all (l2006i) . They 
demonstrated that cloud opacity, adjusted by a sedimentation 
efficiency factor /sed in these models, affects the spectra of all 
dwarfs up to early T types, and the observed CO/CH4 and 
N2/NH3 abundances are indicative of mixing effects equivalent 
to eddy diffusion coefficients between 10^ and 10^ cm^ s"' in all 
atmospheres. But this analysis was still not unable to a unique re- 
lation between spectral type or Tf^g and sedimentation efficiency. 
They also found that the atmospheric parameters derived from 
best fits to individual spectral regions would frequently infer 
different results, or be in disagreement with expectations from 
structural and evolution models. Thus none of the classical static 
models have reproduced the M-L-T spectral transition satisfac- 
torily. 

Attempts have been made to account for atmospheric dy- 
namics in planetary atmospheres, whose models however cat 
not describe the convection cells and the resulting gravity waves 
dMarlev et all 120071: iFortnev et ani2006l) . However, local radia- 
tion hydrodynamics (RHD) models of the surface layers of the 
solar convection have been very successful i n reproducing an d 
analyzing the properties of the granulation (Nordlund ni982h . 
In the meantime, various groups have developed similar codes 
to investigate the atmospheric flows on the sun and other stars 
dSteffen et aV^19^9f: 'Asplund e t al.ll20qo|: ISkartlie n et al.ll2000t 



Stein & Nordlund 2000; Gadu n et al.l |2000; Ro binson etaP 



20031: lVogleni2004) . Amongst others, these models can describe 



self-consistently the mixing of material beyond the classical 
boundaries of a convection zone, as demonstrated fo r instance 
for main-sequence A-type stars jF^vtag et all 1 1996 1) or for M 
dwarfs ^udwig et al. 2002, 20oC 

A treatment of dust within a 
3D s imulatio n of the en velope of an AGB star was included by 
iRevtag & H5fiie^ (I2OO&) . 

The aim of the current work is to extend the latter simula- 
tions into the regime of brown dwarfs, where dust clouds have a 
strong influence on the photospheric temperature structure, and 
to quantify the overshoot from the surface convection zone into 
the atmosphere. 

2. Simulations with C05B0LD 

2.1. Numerical radiation liydrodynamics 

We computed a sequence of 2D RHD models for a gravity of 
10^ cms"^ (loggr=5) and a range of effective temperatures from 
900 K to 2800 K. The models have about 400x300 grid points 
(see Table [1] for details). Most of them are restricted to two 
dimensions because we are unable to cover the prohibitively 
long sedimentation and mixing timescales in 3D: a 2D simu- 
lation in itself takes about one to three CPU-months to com- 
plete. However, a shorter run covering only several dynamical 
timescales and not trying to cover the longer mixing timescales 
is feasible in 3D (mtl5g50mm00n06). 

For this purpose, we used the multi-D RHD code 
CO5BOLE0 (iFrevtag et alJl2002t IWedemever etaHIIOOl in its 



C05B0LD User Manual: http://www.astro.uu.se/~bf/co5bold_main 



local box setup to calculate time-dependent atmosphere models, 
including the very top layers of the convection zone. To realize 
this project, we implemented a dust model (see below) as well 
as dust and low temperature gas opacities. 

The code solves the coupled equations of compressible hy- 
drodynamics and non-local radiation transport on a Cartesian 
grid with a time-explicit scheme. The tabulated equation of state 
accounts for the ionization of hydrogen and helium, and the 
formation of H2 molecules. The ID hydrodynamics fluxes are 
computed with an approximate Riemann solver of Roe-type. 
Because the conditions in the cool objects are almost incom- 
pressible, the fluxes are combined non-split, i.e., the fluxes in 
both the vertical and horizontal directions are computed from 
the same state (and not after each other) and their contributions 
are added. In this way, the generation of spurious pressure waves 
is avoided, which may be produced by a split scheme in regions 
with large gradients but small divergence in the mass flux. 

2.1.1. Dust model 

To account for the presence of dust particles, we added terms 
in the modules for hydrodynamics, radiation transport, source 
terms, and in handling of boundary conditions. It is impossi- 
ble to account for all mi crophysical processes that might play 
a role in dust formation ('Hellin g et alj|200 r; Woitke & Helling! 
I 2OO3I) in current time-dependent multi-dimensional simulations. 
We instead chose a treatment of dust that includes only the most 
important physical processes. The scheme is based o n a simpli- 
fied version of the dust model used in Ho fner et al.l (2003). We 
use a single density field to describe the mass density of dust par- 
ticles and qno5h34 for the monomers (gas constitu ents), instead 
of four for the dust and none for the m onomers as in lHofner et al.l 
( I2003h and'Frevta g & Hofneii ( l2008h . Therefore, the ratio of flie 
sum of dust and monomer densities to the gas dens ity is allowed 
to cha nge, in contrast to the dust description by iHofner et al.l 
(l2003l) . Instead of modeling the nucleation and the detailed evo- 
lution of the number of grains, we assume a constant ratio of 
seeds (dust nuclei) to total number of monomers (in grains or 
free) per cell. If all the material in a grid cell were to be con- 
densed into dust, the grains would have the maximum radius 
'"d.max, which wc havc set to a typical value of 1 jim. This is close 
to the typical particle sizes found for the optically thick part of 
the cloud deck in solar- metallicity brown dwarfs according to the 
DRIFT-PHOENIX models of Witte et al. (2009) and according to 
our own PHOENIX BT-Settl calculations. In both models, parti- 
cle sizes are determined by a balance between settling speed and 
turbulent upmixing according to a basic convective overshoot 
model (cf. iHelling et al .' 2008b for a comparison), and are thus 
in general height-dependent, reaching up to several yum in the 
deepest cloud layers. For the present models, the value chosen 
here should allow a reasonable estimate of the dust opacity in 
the denser parts of the cloud deck. The radius of dust grains 
for given dust mass density pd and monomer mass density is 
computed from 

'"d = '"d.max [Pd/(Pd +Pm)]'^"' ■ (1) 

Condensation and evaporation are modeled as in iHofner et al.l 
(l2003h . parameters and saturation vapor curve adapted to 
forsterite. 

In the hydrodynamics module, monomers and dust densities 
are advected with the gas density. However, according to the 
terminal v elocities given b y the low-Reynolds-number case of 
Eq. (19) in iRossowi d 19781) . a settling speed is added to the ver- 
hhtiy al advection velocity of dust grains, assuming instantaneous 
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model 


P 


n,- X n- 


XXZ 


to 


-t, 


^Cour 


■'in 


rcff.sta 












kmxkm 




1 n't, 

10 s 




-W7- 1 1 

ergK 'g 


K 


K 


K 


mt09g50mm00nl0 


d 


400x430 


220x 94 


90- 


130 


0.30 


0.671 10' 


900 


897 


1000 


mt09g50mm00nl 1 


d 


400x410 


220 x 90 


240 - 


275 


0.30 


0.671 10' 


900 


897 


1000 


mtlOg50mmOOn03 


d 


400x400 


240 X 96 


140 - 


195 


0.27 


0.683 10' 


1000 


1030 


1000 


mtllgSOmmOOnOS 


d 


400x356 


260x 92 


245 - 


285 


0.30 


0.690 10' 


1100 


1114 


1800 


mt 1 2g50mm00n0 1 


s 


400x336 


280x 94 


110 - 


190 


0.40 


0.701 10' 


1200 


1226 


1800 


mtl2g50mm00n07 


s 


400x380 


280x106 


80 - 


150 


0.40 


0.701 10' 


1200 


1228 


1800 


mtl2g50mm00nl0 


s 


400x336 


280x 94 


270 - 


370 


0.37 


0.701 10' 


1200 


1224 


1800 


mtl3g50mm00n01 


s 


400x270 


300x100 


80 - 


190 


0.40 


0.711 10' 


1300 


1335 


1800 


mtl3g50mm00n02 


d 


400x270 


300x100 


130 - 


200 


0.40 


0.711 10' 


1300 


1336 


1800 


mtl3g50mm00n03 


s 


400x270 


300x100 


250 - 


390 


0.35 


0.711 10' 


1300 


1333 


1800 


mtl4g50mm00n01 


s 


400x270 


320x106 


40 - 


140 


0.40 


0.719 10' 


1400 


1436 


1800 


mtl4g50mm00n02 


s 


400x270 


320x106 


140 - 


280 


0.40 


0.719 10' 


1400 


1437 


1800 


mtl5g50mm00n01 


s 


400x270 


340x113 


50- 


140 


0.40 


0.728 10' 


1500 


1533 


1800 


mtl5g50mm00n02 


s 


400x270 


340x113 


150 - 


190 


0.40 


0.728 10' 


1500 


1533 


1800 


mtl5g50mm00n03 


s 


400x270 


340x113 


260 - 


380 


0.30 


0.728 10' 


1500 


1533 


1800 


mtl5g50mm00n04 


s 


300x270 


340x113 


40 - 


150 


0.30 


0.728 10' 


1500 


1532 


1800 


mt 1 5g50mm00n06 


s 


300-X270 


340^x113 


10 


- 15 


0.30 


0.728 10' 


1500 


1532 


1800 


mt 1 5g50mm00n07 


s 


400x300 


340x127 


70 - 


150 


0.30 


0.728 10' 


1500 


1533 


1800 


mtl6g50mm00n06 


s 


400x366 


352x128 


20 - 


190 


0.40 


0.737 10' 


1600 


1648 


1800 


mtl7g50mm00n02 


s 


400x352 


380x133 


40 - 


190 


0.40 


0.747 10' 


1700 


1757 


1800 


mt 1 8g50mm00n07 


s 


400x343 


400x137 


40 - 


190 


0.40 


0.756 10' 


1800 


1858 


1800 


mtl9g50mm00n02 


s 


400x334 


420x140 


40 


- 80 


0.40 


0.765 10' 


1900 


1953 


1800 


mt 1 9g J UmmUUnUi 


s 


400x334 


420X140 


120 - 


290 


0.30 


0.765 10 


1900 


1953 


1800 


mt20g50mm00n05 


s 


400x362 


380x137 


40- 


210 


0.20 


0.775 10' 


2000 


2052 


1800 


mt22g50mm00n05 


s 


400x351 


400x140 


70- 


150 


0.20 


0.798 10' 


2200 


2247 


1800 


mt24g50mm00n01 


s 


400x344 


420x144 


70- 


190 


0.20 


0.825 10' 


2400 


2426 


1800 


mt26g50mm00n01 


s 


400x353 


420x148 


70- 


140 


0.20 


0.859 10' 


2600 


2597 


1800 


mt26g50mm00n02 


s 


400x353 


420x148 


70- 


140 


0.20 


0.859 10' 


2600 


2611 


2800 


mt28g50mm00n01 


s 


400x357 


440x157 


70- 


140 


0.20 


0.905 10' 


2800 


2801 


2800 


mt28g50mm00n02 


s 


400x380 


440x167 


40- 


100 


0.20 


0.905 10' 


2800 


2801 


2800 



The columns show model name, numerical precision (single or double), horizontal x vertical resolution, horizontal x vertical size [kmxkm], time 
span used for averaging [lO^^s], Courant number, entropy of the material in the deeper layers [ergK^'g"'], effective temperature of the PHOENIX 
model used for the start file [K], effective temperature of the RHD model versus the average time span [K], and effective temperature of the 
reference atmosphere used for the opacity table [K]. 



equilibrium between gravitational and viscous forces that act 
onto the grains. 

One problem with modeling the dynamics of dust clouds is 
the span in timescales (short for dust formation and the wave 
period, long for dust settling and thermal relaxation) and spatial 
scales (small-scale dust clouds and possible global flows caused 
by rapid rotation). This is quite similar to simulations of weather 
patterns on Earth, where global wind systems and local cloud 
formation interact. 

Another problem is the poorly known complex micro- 
physics: a complicated chemical network of molecules with 
space- and time-dependent abundances can form dust by means 
of various processes, producing grains with different structures. 
The dynamical behavior and optical properties both depend on 
flie grain type. Furthermore, depletion leads to a change in the 
gas composition that affects the equation of state and gas opaci- 
ties. The current dust model in C05B0LD is designed to repro- 
duce the essential processes, but cannot account for all details 
that might possibly play a role. 

2.1 .2. Equation of state and opacities 

The equation of state accounts for the ionization of hydrogen and 
helium, and the formation of molecular hydrogen. C05B0LD 
can deal with the effects of ionization but not with an ele- 



ment composition that depends on space and time. Therefore, 
the depletion of elements is ignored for the equation of state: 
the formation of molecules has only a minor effect on e.g., the 
heat capacity as long as hydrogen exists in the form of H2. 
However, molecules play a major role for the opacity, and the 
formation of molecules depends both on the abundance and 
depletion of elements. To take this into account, we derived 
the C05B0LD gas phase opacities from monochromatic opac- 
ity tables, k{T, P, v), generated from detailed radiation transfer 
calculations with the general stellar atmosphere code PHOENIX 
( Hauschil dt et al, .1997). We assume full sedimentation of dust 
from the gas phase: the removal of condensable material from 
the gas phase is considered assuming a solar elemental compo- 
sition in f ull phase eqtiilibrium at each tempera ture and pressure 
point (see lAllard et al.l2 001: Ferguson et al. 2005). This approx- 
imation is close to the conditions prevailing in: i) the lower atmo- 
spheric layers that are too hot for dust condensation, ii) the up- 
permost layers where the gravitational settling depletion is par- 
tially compensated by dynamical upwelling of monomers, and 
iii) in the cloud-forming layers as confirmed by observations as 
stated above. The monochromatic gas opacity table was aver- 
aged into 5 bins to minimize the computing time but retain the 
radiative equilibrium properties of the gas. 

In contrast to the sophisticated treatment of the gas opacities, 
we use a simple formula for the dust opacities, which assumes 
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that the large particle limit is valid for all grain sizes and treats 
scattering as true absorption. The dust opacity [cm"'] is 



mt18g50mm00n07 t=1 80060.1 s 



Kd = 3/(2pd,material'"d,max) [PdO^d + Pm)] 



1/3 



(2) 



computed dynamically from the quantities as described for 
Eq. ([T]i in each cell of the simulated atmosphere and added to 
the gas opacity. We concentrate on forsterite grains (Mg2Si04, 
3.3g/cm^) that are relatively abundant and provide the greatest 
contribution to the total dust opacities. 

2.1.3. Boundary conditions 

The side boundaries of the computational domain are peri- 
odic. Usually, an open top is used together with an open bot- 
tom boundary conditions for local models that comprise part 
of a deep convection zone. However, closed boundaries keep 
the amount of dust and condensable material constant within 
the computational domain. We therefore used closed boundaries 
(top and bottom) for all brown dwarf models, although the stel- 
lar convection zone should extend to the center of the star. To 
keep the entropy close to a prescribed value, the internal energy 
is adjusted for a few grid layers (10 km height) at the bottom of 
the model. This mechanism acts as an energy source and replen- 
ishes the radiative energy losses through the top of the model. 
This parameter (the value of the entropy plateau in the deep 
convective layers) controls the effective temperature and is taken 
from the start-up Dusty models. Moreover, a drag force dampens 
downdrafts in these bottom layers. 

The top boundary is closed as well, partly to keep material 
inside. It has a damping zone of about 8 grid points where a 
strong drag force is applied. Damping at an open boundary did 
not appear sufficient to keep gravity waves with moderate Mach 
number (with peak values close to 1) from achieving additional 
growth to implausibly large amphtudes. 

2.1 .4. Initial conditions 

The thermal structure of a start model is based on a classical 
ID stationary stellar atmosphere model produced with PHOENIX 
assuming hydrostatic equilibrium and rad iative plus conv ective 
(using the Mixing-Length Theory, Boh m-Vitensd Il958h flux 
equilibrium. We preferred the dust-rich Dusty over the dust-free 
Cond models even for lower temperatures where the Cond mod- 
els represent dust-free photospheres, because the resulting effec- 
tive temperature of the C05B0LD models agrees very well with 
the effective temperature of the Dusty models (C05B0LD and 
PHOENIX models have per construction the same entropy in the 
deeper layers - not necessarily the same effective temperature). 
We interpolated the Dusty grid points to a finer grid with small 
or no variation in the grid spacing. To allow sufficient volume 
for the surface granules to form, we added several points at the 
bottom by integrating a hydrostatic stratification with constant 
entropy, taken from the bottom point of the Dusty model. At the 
top, we attached a few points with the internal energy value of 
the top point in the Dusty model, to maintain a sufficient distance 
between the top of the cloud layers and the top boundary of the 
computational box. 

We enlarged the model in one horizontal dimension to 
400 points, and imposed small random velocity fluctuations as 
seeds for convective instability. Initially, we set a constant frac- 
tion of the monomers plus dust mass density divided by the gas 
density, but which we reduced somewhat empirically in the up- 
permost layers to account for the partial depletion of material. 




mt18g50mm00n07 t=1 80480.1 s 




As [erg/K/g] 

Fig. 2. Three snapshots of the entropy fluctuations (entropy 
with horizontal average removed) of a brown dwarf model 
mtl8g50mm00n07 with Teff^lSOOK and log g=5. 



The relative amount of material in the monomer bin is then 
determined by the saturation pressure of forsterite. Although 
the Dusty models assume hydrostatic equilibrium, there are 
small deviations from exact numerical equilibrium in the ini- 
tial C05B0LD models. These cause unwanted plane-parallel 
oscillations that we suppressed by a drag force in the initial 
phase of each simulation. To dampen these, we applied a strong 
drag force acting only on plane-parallel motions within the first 
100 sec. In the following 9900 sec, we reduced the drag force 
to remove remaining plane-parallel residuals. For the remainder 
of the run (including the interval where we take averages from), 
we still have a very small but non-zero drag force that dampens 
plane-parallel vertical and horizontal motions on a timescale of 
15 000 sec to suppress some modes that grew in early models 
over very long timescales. 



3. Results of the simulations 

3.1. Convection and timescales 

The quiet solar surface - far away from sun spots - is character- 
ized by a mottled pattern of bright hot rising areas surrounded 
by dark lanes of cool downflowing material - the so-called gran- 
ulation at the top of the solar convection zone. Because of their 
higher surface gravity and lower effective temperature (log g=5, 
Teff^lSOOK), granules on brown dwarfs are about a factor 
of 10 smaller than their counterparts on the sun (logg(=4.44, 
T^eif =5775 K). Only tiny velocities are required to transport the 
energy flux through the convection zone, resulting in nearly in- 



Freytag et al.: Convection, overshoot, and gravity waves in low-mass dwarfs 



5 




100 200 300 

X [km] 



Fig. 1. This snapshot from a brown dwarf simulation with reff=1858K, log g-5 shows the velocity field as pseudo-streamlines, 
color-coded according to the dust concentration. The flow in the lower part is due to the surface granulation of the stellar convection 
zone. The top is dominated by gravity waves. 




X [km] 

Fig. 5. As Fig. [Jfor a brown dwarf model (mtlOg50mmOOn03) with Teff^lOOO K and log g^5. 



compressible low-Mach-number flows (with a maximum verti- 
cal rms value ranging from about 0.06 at reff=2800K to about 
0.003 at reff=900K). Accordingly, only small-amplitude pres- 
sure waves are present in the brown dwarf models. For compar- 
ison, the acoustic timescale (for an up-down-up wave travel) is 
about 90 s for a 2800 K model and about 60 s at the cool end of 
our sequence. The Brunt-Vaisala period in the stable layers is 
about 25 s, while the convective growth time in the unstable lay- 
ers increases with decreasing eff'ective temperature from around 
40 s to about 300 s (minimum values at the "most unstable lay- 
ers"). The free-fall timescale (to drop one pressure scale height) 
is around 3 s in the photosphere. 

To check the transition from the starting conditions to a 
quasi- stationary state, we consider time sequences of spatial- 
averaged quantities, such as temperature, rms velocities, and 
dust concentration. Starting from random initial fluctuations, the 
onset of convection takes a few 100 seconds, longer at lower 
effective temperatures. A statistically stable pattern develops af- 
ter a few 1000 seconds. Wave amplitudes relax on somewhat 
longer timescales. In contrast, the thermal relaxation time - 
particularly of the deeper convective layers - is much longer. 
However, the thermal structure of these layers - an adiabat - 
remains essentially the same as in the initial model due to our 



choice of treatment of the lower boundary (keeping the entropy 
constant instead of imposing a certain flux). In this way, there 
is no need to cover the complete thermal relaxation time. The 
longest timescale to be covered is the relaxation time of the dust 
concentration that has an effect onto the temperature structure. 
Therefore, each simulation covers a few days of stellar time. The 
hydrodynamic time step is about 0.03 s and because of the rela- 
tively long radiative relaxation time, we perform multiple (typi- 
cally 6) hydrodynamical sub-steps per radiation transport step. 

Snapshots from our atmosphere simulations are presented in 
Figs.[T]through[8]while the complete videos are provided as sup- 
porting materiaQ. Figures[T]and|5]use pseudo-streamlines to vi- 
sualize the flow field. Figures |2] (for a 1800K model) and|6](for 
a 1000 K model) display sequences of the typical granulation 
pattern, cool downdrafts occurring in a warmer environment. It 
is clearly separated from the atmosphere in the upper half of 
the box that shows inhomogeneities induced by gravity waves. 
The downdrafts are relatively narrower than in solar granulation 
(Ludwig et al. 2002, 2006). In the image sequences, the first pair 
is 20 s apart whereas the last snapshot is taken several minutes 
later 



^ |http://phoenix.ens-lyon.fr/papers/FreytagEtA12009/| 
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Fig. 3. Three snapshots of the dust concentration of a brown 
dwarf model mtl8g50mm00n07 with refF=1800K and log g=5. 



Fig. 4. Three snapshots of the concentration of dust-fmonomers 
of a brown dwarf model mtl8g50mm00n07 with reff=1800K 
and logg(=5. 



The entropy profiles - averaged horizontally over constant 
height and in time — in Fig. |9] (top left) show a strong increase 
in the upper atmosphere, with only a minor drop at the top of 
the convection zone, and an almost flat distribution inside the 
convection zone. This is indicative of very eflicient convection 
resembling typical conditions in the stellar interior. 

3.2. Exponential overshoot 

The typical magnitude of the velocity fields can be inferred from 
the plot of the rms of the vertical velocity versus pressure for var- 
ious eff'ective temperatures in Fig. |9] (middle panels). The con- 
vective velocities fall significantly from the peak value inside the 
convection zone (on the right) until the top of the unstable lay- 
ers, and even further into the overshoot region. Th e scale height 
of ex ponenti ally decreasin g overshoot velocities (iFrevtag et al.l 
[199 6: Ludwi getal.ll2006h is so small that they do not induce 
significant mixing in the cloud layers about two pressure scale 
heights further up. Nevertheless, they are able to mix material 
across the boundary between stable and unstable layers. 

The bottom right panel in Fig. [TO] shows the relatively large 
scale height of the convective velocity at high effective tempera- 
tures. This extended overshoot may play a role in the replenish- 
ment of dust material. However, the overshoot scale height de- 
creases rapidly with Tgff and remains small (//„»0.28 Hp) from 
about 2200 K on, indicating that this type of overshoot is in- 
significant for material mixing within the forsterite cloud lay- 
ers. On the other hand, for dust types that form at s lightly higher 
temperatures (around 2000 K) as discussed e.g., in lHeUing et al.l 



(l2004h the mixing caused by convective overshoot can play a 
role. 

3.3. Gravity waves 

It is instead gravity waves that dominate the mixing of the at- 
mospheric layers (the upper half of the models in Figs. [T]to [8]) 
with periods of about 30 to 100 seconds and amplitudes that in- 
crease with height (Fig.|9]l. Most prominent is the fundamental 
g-mode, visible particularly in Fig.|6]as a significant brightening 
between heights of and 20 km. In addition, there are several 
modes with larger horizontal and vertical wave number. These 
waves show up together with the first surface granules, well be- 
fore the downdrafts "hit" the lower boundary. This indicates that 
the granular flow as such is responsible for the wave excitation, 
and not artifacts related to the way flows at the lower boundary 
are handled. 

Figures [T] and [3] demonstrate the location of the dust clouds 
and the effect of the thermal inhomogeneities induced by the 
gravity waves onto the dust concentration. The generated small 
amount of vertical mixing (the wave motion is mostly reversible) 
is sufficient to balance gravitational settling of dust grains and 
allow dust clouds to form in the hotter models. In addition, dust 
concentration and cloud thickness are modulated by the waves 
because of the induced temperature fluctuations. 

Atmospheric gravity waves are a common phenomenon. On 
Earth, they are known to form clouds over e.g., the US mid- 
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Fig. 6. Three snapshots of the entropy fluctuations of a brown Fig. 7. Three snapshots of the dust concentration of a brown 
dwarf model mtlOg50mmOOn03 with reif=1000K and log g=5. dwarf model mtlOg50mmOOn03 wifli Teff^lOOOK and logg=5. 



west plainfl Their energy release is involved in heating the ex- 
ospheres of Jovian planets as observed for Jupiter from Galileo 
probe results (Young 1998). 

Simulations of convection producing g ravity wav es in stellar 
conditions have a long tradition (Hurlburt et al.lll98a) . However, 
the quantitative estimate of the amplitude and the true detec- 
tion of internal gravity wav es can be a difficult ta sk, even for 
the well-studied solar case jBelkacem et al.ll2009l) . The detec- 
tio n of gravit y modes that probe the solar core was announced 
bv lGarcia et a l. (2007). 

The initial phases of simulation indicate that gravity waves 
are generated near the top of the convection zone (see e.g., 
iDintrans et al.ll2005l) . The gravity waves are produced by non- 
stationary downdrafts "sucking" at the stable photospheric lay- 
ers. In this way, the downdrafts are able to inject kinetic en- 
ergy into the photosphere and to transport some material from 
there into the deeper convection zone. However, there are no 
obvious "events" o f wave generation as for p -modes in the sun 
(iGoode etal.lll998t IStein & Nordlund"200D or in the simula- 
tions of gravity waves generated by an idealized co nvection zone 
embedded between stable layers by IDintrans et al.l (12005). 

The mixing efliciency of the waves increases rapidly with 
height - steeper than expected from the mere growths in ampli- 
tude caused by the increasing non-linearity. This could be the 
dynamical updraft mechanism responsible for the upwelling of 
N2 and CO gas observed via the enrichment of CO and deple- 



Mesonet, 



I. 



E. 



2007, 

|http://mesonet.agron.iastate.edu/cool/l 



Gravity Wave Movie: 



tion of CH4 and NH3 absorption bands in the spectra of T dwar fs 
dSaumon et al ] |2006t [Stephens et alJl2009HGeballe et alll2009h . 



3.4. Convection within dust clouds 

The fluctuations in the dust concentration in the 1800K model 
in Fig. [1] are mainly induced by up and down motions of gravity 
waves that provide an inefficient mixing that balances the set- 
tling of dust grains. However, when the optical thickness of the 
dust clouds becomes sufficiently high, convective motions within 
the dust clouds start to develop and provide more efficient mix- 
ing of material (cf. the dust concentration of the lOOOK model 
in Fig.|7]i. However, in the snapshots the fluctuations and flows 
due to the waves somewhat obscure the dust cloud convection, 
whereas the overturning motions are clearly visible in movies 
and have a different signal in a k-u) diagram. 

There are different intermittent processes: occasionally, ma- 
terial from the dust layers is dredged up to the layers with rel- 
atively low dust concentrations above the clouds. The grains 
quickly fall back. But monomers can remain a while, until they 
condense into dust at the top of the cloud deck. The cloud layer 
thickness varies not only with the wave on a timescale below 
one minute but also in irregular cycles on timescales of hours. 
The irregularity and amplitude increases with decreasing effec- 
tive temperature. 

During the initial phases of a simulations, a violent thin 
cloud convection zone develops for a limited time until the 
model is relaxed. This phenomenon relates to differences be- 
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and log g-5. 



tween our start model and the final outcome. However, on ac- 
tual brown dwarfs large-scale flows might cause an imbalance 
in the local dust concentration that leads to a similar localized 
enhanced cloud activity. 

3.5. Dust and stratification 

In the top left panel in Fig[TOl we show the location of the dust 
clouds (circles connected by vertical lines) relative to the under- 
lying gas convection zone (located below the crosses). 

The mixing processes within the dust cloud layers have dif- 
ferent height regimes. At the bottom of the dust clouds, the tem- 
perature varies around the condensation value but there is lit- 
tle mixing. With our dust scheme, which assumes the presence 
of nuclei everywhere where dust or monomers are present, dust 
forms and evaporates during these temperature cycles (see the 
dust concentration at z~10km in Fig. |7]). The dust formation 
would be more difficult if new dust grains had to nucleate, be- 
cause that would require some level of supersaturation. Within 
the clouds, material is mixed by gravity waves and/or convec- 
tion (depending on eff'ective temperature). The top of the clouds 
is sharp but inhomogeneous due to (sometimes braking) waves 
and cloud convection. Above the cloud layers, there are still mix- 
ing flows that try to equalize the concentration of monomers 
with height. The concentration value depends on the efficiency 
of mixing, dust formation, and dust settling in the cloud layers 
below (Fig.|9] bottom panels). 



Dust clouds have a strong effect on the thermal structure 
(Fig. |9] top right panel): there is a fairly shallow temperature 
slope beneath the cloud layers with values of about 1600K be- 
cause of the greenhouse effect, a rapid drop within the clouds due 
to the large dust opacities, that can even drive cloud convection, 
low temperatures (with values around lOOOK and small varia- 
tions) in the mostly dust-free upper atmosphere, and in some 
cases a small increase at the top of the models of about 100 K be- 
cause of the dissipation of kinetic wave energy. At some height 
above the cloud, gravitational settling of dust grains becomes 
more efficient than mixing. The dust density drops rapidly and 
with it dust opacity and temperature, causing a rather sharp (but 
variable in space and time) upper boundary of the clouds. The 
concentration of dust and monomers (material that potentially 
can form dust) in Figs. |4] and [8] shows complete mixing in the 
convection zone, depleted layers at the top of the atmosphere 
(due to gravitational settling), and a partially mixed region in- 
between. 
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Fig. 9. Various averaged quantities versus logarithm of pressure: Top left: mean entropy for various effective temperatures and 
loggr=5. The plus signs mark the layers with Rosseland optical depth 10"^. The legend is the same in all panels. Top right: mean 
temperature. Middle left: logarithm of rms horizontal velocity. Middle right: logarithm of rms vertical velocity. Bottom left: loga- 
rithm of monomer concentration. Bottom right: logarithm of dust concentration. 
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3.6. Effective temperature dependency 

We summarize the dependence of our model properties on effec- 
tive temperature in Fig. [10] The upper right panel shows that 
the atmospheric velocities (squares) do not follow the mono- 
tonic decrease in the convective velocities (crosses) up to the 
lowest effective temperatures, but rise after a minimum around 
^eff = 1900 K. However, the cloud thickness (top left panel) and 
cloud mass (bottom left panel) increase monotonically with de- 
creasing effective temperature. The cloud extension with Teff is 
slightly erratic because the models have not perfectly converged. 
At higher effective temperatures, the thin see-through clouds al- 
low a view of the upper layers of the convection zone (plus signs 
in both top panels), while at the lower- Teff end most of the cloud 
mass sits below the visible layers. 

The rms vertical velocity, which is often used in static model 
atmospheres to estimate line broadening, non-equilibrium chem- 
istry due to mixing, or cloud formation, can be parametrized 
as a function of effective temperature for our sequence of 30 
models. Hence, with the logarithmic normalized temperature 
X = log(min(max(r, 900 K), 2800 K)/l^:), we obtain for the ve- 
locity scale height of the wave amplitude (bottom right panel in 
Fig. [loll 

HJHp = max(2.855, -43.1 4- 14.34 x) , (3) 

and for the logarithmic ratio of maximum convective velocity to 
wave amplitude extrapolated to this layer 

logr„ = max(835.684 - 828.156;c + 273.426x2 - 30.1148x^ 

-25.496 + 7.1 104x) , (4) 

where the first expression in the max functions is a good fit for 
temperatures below approximately 2000 K, and the second ex- 
pression for temperatures above. To recover a crude estimation 
of the mixing efficiency based on these formulae, we propose to 
compute the amplitude Vmax and position Pmax of the maximum 
convective velocity with e.g., the Mixing-Length Theory, to add 
log r„ to the logarithmic velocity amplitude, and to extrapolate 
from this starting point the wave amplitude with H^/Hp into the 
atmosphere, 

log V = log V^,, + log r„ - (log P - log P^,.) I {HJHp) , (5) 

where P is the pressure. We note, that for the models with 
reff=900K and 1000 K the velocities drop after reaching a max- 
imum of about 10ms while the rise in the velocities of the 
hotter models appears to be limited only by the top of the com- 
putational domain. 

The rise of velocities above the convection zone becomes 
steeper with decreasing temperature (bottom right panel in 
Fig.[TOb. But the amplitudes in the atmosphere (Figs. l9l and [TOl 
top right) fall from models with high effective temperatures to 
models with T^s ^ 2000 K. 

For even cooler models, the amplitudes increase again. This 
increase is not because of the velocities in the underlying gas 
convection zone that decline steadily as effective temperature de- 
creases. Instead, below 2000 K, the clouds have grown to such a 
large vertical thickness and density that cloud convection begins 
- with effects onto the atmospheric velocities and temperature 
structure that increase with further decreasing effective tempera- 
ture. We attribute the rise in velocity for low-temperature models 
mainly to the emergence and growth of cloud convection. 

Therefore, the dust clouds also affect the waves: when they 
become extended enough they split the atmospheric cavity for 
gravity waves into two separate zones by generating an entropy 



plateau within the atmosphere (Fig. |9]l- At T^s below about 
lOOOK, the gravity waves are trapped mainly inside the region 
between the two convection zones. The layers above the dust 
convection zone exhibit only small wave amplitudes. 

The convective-radiative boundary becomes steeper, hence 
harder, with lower effective temperature, easing the gravity wave 
generation. In addition, there is a slight change in the topology 
of granules: in the hotter models, the downdrafts that delimit 
the granules are of roughly similar strength and merge occa- 
sionally, while in the cooler models just a few (2 or 3) "super 
downdrafts" dominate and absorb the smaller ones that form on 
top of the granules. This process occurs with higher frequency 
than would be expected if the merging type were the same as in 
the hot models - with possible consequences for the interaction 
between convection and waves. 



4. Discussion 

4. 1 . Parameter dependence 

The details of the wave generation process and the type, am- 
plitude, interaction, long-term evolution, and spectrum of the 
waves are complex and may have connections to numerical set- 
tings. Some small-amplitude pressure waves in the hottest mod- 
els are emitted from the non-stationary downdrafts as expected 
and become invisible at intermediate temperatures. However, in 
the lowest-temperature models (reff<1300K) the signature of 
high-frequency p-modes unexpectedly showed up in the pressure 
fluctuations, where g-modes and convection have only a small 
signal due to their nearly incompressible nature and small am- 
plitude. The p-modes contribute to the velocities only close to 
the very top of the computational box. Their amplitude depends 
sensitively on viscosity and the position of the top boundary. 
They vanish, when the Courant number is reduced from 0.4 to 
about 0.3, depending on other model details (the usual stability 
criterion sets an upper limit at 0.5 for the 2D models). The mod- 
els that we used for our final analysis show no or only traces of 
these p-modes. 

An early version of the models showed (in addition to the 
"normal" spectrum of gravity waves that occur as soon as con- 
vection sets in) after the simulation had run for a long time a 
slowly exponentially growing gravity wave in the fundamental 
mode. It grew until the code crashed because of too steep veloc- 
ity gradients at the top of the box. It had relatively little effect 
on mixing, but induced temperature fluctuations modulating the 
dust concentration. Limiting the model depth, and using both 
a finer vertical grid and a smaller Courant number prevented an 
exponential growth of the mode. However, the mode itself is still 
present and quite prominent in the cooler models. 

We varied several numerical parameters to check their in- 
fluence. By considering a pair of models with reif=1300K, one 
in single, the other in double precision, we found hardly any 
difference at all in the resulting mean properties (velocities and 
temperature). However, at reff<1100K the density fluctuations 
in the convection zone become so small that they cannot be re- 
solved properly using single precision: after a transient phase, 
convection dies out leaving only small-scale low-amplitude ve- 
locity fluctuations in the "convection zone" that are due to round- 
off errors. Therefore, all runs with reff<l lOOK were performed 
in double precision from the initial tests on. 

For a 1500 K model, we decreased our standard horizontal 
resolution by going from 400 to 300 horizontal grid points and 
found no noticeable difference in the mean structures, although 
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Fig. 11. Mean temperature versus logarithm of pressure for two 
models with Teft ~2600 K, log g=5, but different opacity tables. 



ment at reif=900K, 1 200 K and, 2800 K with similar agreement 
between the pairs of curves. 

Simulations that are not yet complete and that will be pre- 
sented in future publications include a sequence with other 
gravity values that shows no qualitative change in the out- 
come, although convective velocities, wave amplitudes, and 
dust-formation rate equations noticeably depend on gravity: the 
dependence of the flow field and the cloud thickness on eff'ective 
temperature will be different at other gravities. 

A 3D model with reff= 1500 K (after taking about half a year 
to cover 4 hours of stellar time with six CPUs) has completed the 
transition from the initial 2D configuration to a fully 3D flow pat- 
tern. We see only small changes in the wave amplitudes and con- 
vective velocities relative to a 2D run (Fig. [T2]) . The 3D model 
data are in close agreement with the 2D models in the scatter 
plots in Fig (TO] Only its overshoot scale height (bottom right 
panel, Hp) is noticeably smaller than the average of the 

2D models at reff=1500K (//„=0.3 Hp). 

The exploration of other parameters such as grain size and 
other types of dust awaits further simulations, using a more de- 
tailed cloud model (multi-size-bin scheme, in preparation). 
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in the geometry (horizontal grid points, vertical extent, dimen- 
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the thin convective downdrafts and some small-scale cloud struc- 
tures are somewhat less well resolved. 

For the generation of our binned opacity tables, the bins were 
optimized for three reference atmosphere cases (reff=1000K, 
1800K, and 2800K with log0=5.O and solar metalhcity). The 
resulting tables are most accurate for model parameters close 
to these values. In Fig. (TT] we show the temperature struc- 
ture of two models with reff=2600K using an opacity table 
made with a 1 800 K reference atmosphere in mt26g50mm00n0 1 , 
and with a more appropriate 2800 K reference atmosphere in 
mt26g50mm00n02. The resulting differences are most important 
for the outermost layers, but remain negligible for determining 
the mixing mechanisms in these atmospheres. 

Figure [T2] shows the rms velocity for four 1500 K runs with 
different positions of the upper boundary. As expected, the more 
extended model has a larger peak velocity, while the agreement 
between the two curves is good in the lower atmospheric layers 
and excellent in the convection zone. We tried the same experi- 



4.2. Comparison with previous simulations 

iLudwig et al] (l2002l |2006) studied the structure of mid M- 
type, dust-free atmospheres including the mixing properties of 
macroscopic flows with 3D hydrodynamical simulations. Our 
reif~2800K model coincides in temperature with the coolest 
model of Ludwig et al., and this allows a comparison with the 
results. While the maximum rms vertical velocities in the con- 
vectively unstable layers turn out to be similar, the atmospheric- 
wave-dominated velocities are about 50 % higher in Ludwig et 
al. Importantly, the scale height of the decline in the convective 
velocity field amplitude into the stably stratified layers is found 
to be similar (H^/Hp = 1.2 this work, 1 . 1 in the work of Ludwig 
et al.). At first sight, this may appear surprising considering the 
systematic differences expected between simulations conducted 
in 2D and 3D as discussed by Ludwig & Nordlund! (l2000h : in 
2D the efficiency of wave generation is usually higher, and the 
transition between convectively stable and unstable regions is 
more gradual. However, the work of Ludwig & Nordlund refers 
to higher Mach-number flows, and pressure waves, not gravity 
waves, which are relevant here. In the higher Mach-number con- 
ditions studied by Ludwig & Nordlund, towards lower Mach- 
numbers the sharpness in the stable-unstable transition becomes 
more similar in the 2D and 3D simulations, so that the similarity 
in the M dwarf regime appears plausible. This is also borne out 
by a comparison with a 3D model compiled for 2800 K, which 
has a very similar velocity profile to its 2D counterpart. 

Ludwig and collaborators had reasons to believe that the 
gravity waves present in their models were an artifact of the 
lower boundary condition. Moreover, they argued that the mix- 
ing efficiency is too small - because of both low (i.e., linear) 
amplitude and the insufficient shearing - to produce small-scale 
turbulence due to Kelvin-Helmholtz instabilities. Since in their 
models, convective overshoot was potentially able - when ex- 
trapolated to lower Teff - to keep dust grains in the optically thin 
layers, they took assumed that waves are not important and con- 
vective overshooting is sufficient to explain the presence of dust 
clouds in brown dwarfs. 

The present calculations cover the actual parameter regime 
of dust harbouring atmospheres. They show that - in contrast to 
expectations motivated by hotter models - convective overshoot 
alone is not capable of keeping dust grains in the atmosphere. 
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Overshooting motions decline more rapidly towards lower effec- 
tive temperatures (Fig.fTOll. Our thorough investigation of the in- 
fluence of the boundary conditions on the excitation of the grav- 
ity waves indicate that they are indeed intrinsic to the flow evo- 
lution proper and not a numerical artifact. The gravity waves's 
ability to mix is indeed low but the waves remain in our hotter 
local models nevertheless the most efficient process, accompa- 
nied by dust convection in the case of heavy dust formation. All 
in all, we consider our present res ults consistent with the find- 
ings of iLudwig et al.l (l2002i l2006l) . but reassign the importance 
to the mixing by waves. 

A comp lementary approach to ours is pursued by 
iHelling et al.' (2004): rather than on macroscopic scales (pres- 
sure scale heights, depth of the atmosphere, granular diameter), 
they concentrate on mesoscopic scales. Their 2D model is about 
as large (500x500 m^) as one of our grid cells. They investigate 
the influence of driven turbulence, represented by a set of im- 
posed pressure waves, on the formation of dust, particularly in 
regions where the temperature is slightly too high (T=2100K) to 
allow nucleation in an undisturbed atmosphere. We agree with 
their findings that fluctuations in the thermodynamic quantities 
can have an influence on the dust formation process. However, 
we identify gravity waves and not pressure waves as important 
contributors to the mixing in BD and M dwarfs, in addition to 
with convection within thick clouds and convective overshoot 
very close to the underlying gas convection zone. An important 
parameter in their simulations is the Mach number of the induced 
acoustic waves, for which they assume values of about 0.1 in ID 
models and 1 in 2D models. However, peak convective Mach 
numbers (taking vertical and horizontal velocities into account) 
in our models are between 0.1 and 0.01, and rapidly decrease in 
the overshoot regions where high-temperature dust might form. 
The amplitude of turbulent structures on the grid cell scale and 
below - that we obviously cannot resolve in our models - would 
be even smaller. And only a tiny fraction of the energy can be ex- 
pected to be transformed into pressure waves under these nearly 
incompressible low-Mach-number conditions. Therefore, based 
on our simulations we cannot justify the assumption of almost 
sonic p ressure waves in the atmospheres of brown dwarfs as 
made in lHelhng etaP (l2004 . 
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Fig. 13. Logarithm of diffusion coefficient versus logarithm of 
pressure for a model with refF=1000K according to Eqs. (|6]l 
to The black crossed indicated estimates derived from the 
the averaged vertical fluxes and concentrations of dust and 
monomers. 
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4.3. Diffusion coefficient estimate 

One can model the mixing of material by macroscopic flows - 
on average - as a diffusion process. However, in the hydrody- 
namical models there is a correlation between the sign of the 
vertical motions (upward of downward) and the grain growth. 
At the same velocity amplitude, the mixing efficiency of con- 
vective overturning flows is also much higher than that of (nearly 
reversible) wave motions, causing errors in the translation from 
the rms velocities to the actual mixing efficiency. The mixing ef- 
ficiency can however be estimated from the rms vertical velocity 
of our model sequence as in Eq. (|5]l. And the diffusion coeffi- 
cient can be estimated from the local vertical velocity and the 
pressure scale height Hp as typical length scale via 



D cc VH„ 



(6) 



However, the waves have a varying amplitude with height and 
therefore the typical length scale is not constant. On the other 
hand, their period is rather close to the Brunt-Vaisala period, 
which can be used as a characteristic timescale. Using the similar 
acoustic period, = 2//p/csound, we obtain 



log(P/[dyn/cm']) 

Fig. 14. Logarithm of diffusion coefficient D (or K-^) according 
to Eq. ^ versus logarithm of pressure for the same set of models 
as in Fig.|9] 



where the Mach number Ma = V/csound- To take into account 
the increase in mixing with increase in non-linearity, one could 
multiply with the Mach number again to obtain 



DocMa^y//,, 



(8) 



Docy2f ccMaV//,, 



(7) 



Profiles for the diffusion coefficients according to the Eqs. (|6]l — 
dHJ (replacing "oc" by "=") are displayed in Fig. [13] Additional 
crosses mark estimates based on the horizontal and tempo- 
ral averages of vertical flux and density profiles of dust plus 
monomers. Because the flux is divided by the vertical deriva- 
tive of the concentration, which can be very small, these values 
are not well-behaved everywhere. 

The diffusion coefficient in brown dwarfs is not a "constant 
of nature" but depends on the physical process driving the mix- 
ing, the height in the atmosphere, and the effective temperature. 
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Fig. 15. Logarithm of relative bolometric flux intensity contrast 
plotted versus effective temperature for all models in Table [1] 
Top curve (circles): total contrast (spatial plus temporal contri- 
bution). Bottom curve (plus signs); temporal variations only. 



A lower gravity would lead (via increased convective velocities 
and larger time and spatial scales) to high diffusion coeflicients. 

Scaling relations such as Eqs. (|7]l and (O can serve as a 
first step in describing the diff'usion properties in ultracool atmo- 
spheres in greater detail. In addition, they can easily be translated 
into classical atmosphere codes. In an earlier study, we have im- 
plemented a non-equilibrium chemistry model in the PHOENIX 
BT-Settl code, using the diffusion coefficients derived from the 
overshoot contribution of the RHD simulations only. In atmo- 
sphere models of the Tl brown dwarf e Indi Ba based on this 
approach, we obtain characteristic values of 10^ < D < 10^ 
cm^ s"' for the transition region from CO- to CH4-dominated 
carbon chemistry. We found the resulting non-equilibrium abun- 
dances of CO in the line-forming region, using an upda ted and 
more efficient reaction model than lSaumon et al.] (l2006l) . to just 
slightly underestimate the observed CO line strengths in this 
benchmark transition dwarf (King et al. 2009). This indicates 
that somewhat more efficient mixing than provided by over- 
shoot alone is required, thus supporting an additional contribu- 
tion from gravity waves. The latest revision of the BT-Settl mod- 
els beeing tested to include the effect of both dust mixing and CE 
departures, also find that Eq. O provides a close match to ob- 
servational constraints. 



In a more phenomeno logical approach lSaumon et al ] d2006h 
and lStephens et al.](l2009l) explored the effects of a constant eddy 
diffusion coefficient above the Schwarzschild boundary on ob- 
servable departures from nitrogen and carbon equilibrium chem- 
istry. They found best-fit values ranging from 10- to 10'' cm s , 
although based on slower reaction rates (see above). Allowing 
for the differences in the reaction scheme, their results thus agree 
with the range of values that we find for the diffusion coefficient 
(crosses in Fig. [T3l l in the region spanning the top of the over- 
shoot layers, the gravity-wave region, and the base of the cloud 
layers (7.2 > logP > 5.5). Within the clouds, we find larger 
values. 



4.4. Brightness variations 

Our simulations show temporal intensity variations for a wide 
range of timescales from half a minute to hours. We find low- 
amplitude (less than 1%, see Fig.fTSl) and short-period variability 
(1 min) due to the gravity waves producing temperature fluctua- 
tions and modulating the dust density and the vertical thickness 
of the clouds. Relatively short-lived (several minute) phenomena 
are the occasional dredge-up or outburst of material above the 
clouds, where dust falls back rapidly while monomers remain 
much longer. In the coolest models, the gravity waves are less 
visible in the intensity fluctuations, which are dominated instead 
by aperiodic variations on the scale of hours with an amplitude 
of a few per cent. But these results barely reach the scales of the 
observed variabifity of L-type brown dwarfs, which is often ape- 
riodic (scal es of hours to days) a nd of low amplitude (mmags) as 
reported by iGelino et al.l (l2002h . Still, the spatial intensity con- 
trast found in the models (Fig. [TSl l with significant dust layers 
is significantly higher than the contrast that would be induced 
by granulation alone (the small contrast for models above about 
2500 K). 

Our results are only those for a patch of the surface and the 
variability amplitudes obtained will average over the rest of the 
surface. To determine the brown dwarf surface distribution of 
clouds, one must go beyond the present simulations to 3D mod- 
els that are as large as possible and include rotation effects. 

We have neglected the effects of rotation despite the rapid 
rotational periods of brown dwarfs (P < 4 hrs). The convec- 
tive turnover time in the box is several minutes, which is short 
in comparison. Neither the surface granules, nor our rms veloc- 
ities, should therefore be severely affected. But granules could 
move with a global meridional flow. Other global flows caused 
by rapid rotation may exist that could move the dust around. 

A cloud cover disruption has indeed been suggested as 
a possible additional cause - the cloud layers sinking rela- 
tive to the fine formin g layers - of the L-T spectral transition 
dAckerman & M arlev'2001) that could lead to weather phenom- 
ena and spectroscopic variability. 



5. Conclusions 

We have performed radiation hydrodynamics simulations with 
C05B0LD of a sequence of brown dwarf atmospheres extend- 
ing previous studies to lower temperatures. The numerical model 
includes a simple treatment of the formation and destruction of 
dust, as well as its gravitational settling and advection, and also 
the interaction with the radiation field. 

We provide a fit to the rms velocity in the atmosphere that 
can be used to estimate the mixing. The convective velocities 
fall significantly from the peak value inside the convection zone 
to the top of the unstable layers, and even further into the over- 
shooting region. However, the scale height of exponentially de- 
creasing overshoot velocities is so small that they do not induce 
significant mixing in the cloud layers. Above a local minimum in 
the vertical velocities, gravity waves dominate in the hotter mod- 
els with an amplitude and mixing efficiency that increase rapidly 
with height, enough to balance the gravitational settling of dust. 
The wave amplitude decreases with decreasing effective temper- 
ature. In the cooler models, the dust layers are thick enough to 
cause convection within the clouds leading to efficient mixing 
within the cloud layers. 

Models with high effective temperatures (2500 K < T^ff < 
2800 K) show a high-altitude haze of optically thin forsterite 
clouds. At lower effective temperatures (Teff < 1400 K), thick 
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and dense forsterite clouds exist but mostly below the visible 
layers, which are essentially depleted of the material that went 
into the dust. For intermediate effective temperatures, dust is an 
important opac ity source in the atmosp here. This agrees with 
observations bv lGolimowski et al] (|2004|) . which place i) the on- 
set of important refractory element depletion, where both TiO 
and VO bands weaken in spectra because of condensation of ti- 
tanium and vanadium, and greenhouse effects at about 2500 K, 
ii) the maximum greenhouse effects at about 1800K (M to L 
transition), and iii) the transition between dust-rich and dust-free 
brown dwarfs (L to T transition) at around 1450 K. We there- 
fore feel confident that the mixing efficiency determined by our 
simulations is adequate. Although an investigation of the spec- 
tral properties of the models exceeds the scope of this paper, the 
formulae that we provide for the velocity field will allow the dis- 
crimination between diverse cloud model assumptions for brown 
dwarfs and planetary atmospheres. 
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